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ABSTRACT 


A linear triangular finite element formulation was used 
to solve the steady state two dimensional conduction heat 
transfer equation. A FORTRAN IV computer program of the 
above formulation, employing double precision arithmetic 
and compact storage techniques, was applied to study heat 
transfer in an internally finned rotating heat pipe. Results 
were obtained for water in a copper and stainless steel 
condenser with varying outside heat transfer coefficient, 
rotational speed, and fin angle (number of fins). 

Numerical results of the heat transfer rate in the copper 
condenser were shown to have less than 2% variance to those 
obtained earlier. A rotating heat pipe designed to operate 
With a small value of outside heat transfer coefficient 
E Orrences no significant increase of the heat transfer rate 
with an increase in rotational speed, since the value of the 
outside thermal resistance dominates. 

Heat transfer rate continuously increases as the fin 
angle decreases. However, the increase is only slight when 


the fin angle is less than ll degrees. 
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te etl TRODUCTION 


ни THE ROTATING HEAT PIPE 

The rotating heat pipe is a closed container designed 
ЕЕ стапсгег а large amount of heat in rotating machinery. 

Its three main component parts are: a cylindrical evaporator, 
a truncated cone condenser, and a working fluid as shown 
Ші Figure l. 

At rotation above the critical speed of a rotating heat 
pipe, the working fluid forms an annulus in the evaporator, 
and will be vaporized by heat addition in it. The vapor 
flows toward the condenser, as a result of a pressure differ- 
ence, transporting the latent heat of evaporation with it. 
External cooling of the condenser causes the vapor on the 
inner wall to condense and release its latent heat of evapor- 
ation. The centrifugal force due to the rotation has a 
component acting along the condenser wall that will act to 
drive the condensate back to the evaporator where the cycle 
is repeated. 

In a conventional heat pipe, the force driving the con- 
densate back to the evaporator is due to capillary action, 
which poses a limit to its operation. However, the driving 
force in the rotating heat pipe can be designed to operate 


Meany Orientation, including in a zero gravity field. 


159 
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O TILIZATION OF AN INTERNALLY FINNED CONDENSER 

The first theoretical investigation of the rotating 
heat pipe conducted at the Naval Postgraduate School, was 
performed by Ballback [1] in 1969. He studied the limi- 
tations of performance imposed on the rotating heat pipe 
due to various fluid dynamic mechanisms. Using existing 
theory and experimental correlations, he was able to esti- 
mate the sonic limit, boiling limit, entrainment limit, 
and the condensing limit of performance. 

Tne Sonic Limit 

When increasing the heat flux in a rotating heat 

Pipe, it is entirely conceivable to reach a limiting flow 
rate of the vapor brought on by a choked flow condition in 
the pipe. This condition imposes a limiting value on the 
amount of energy the vapor can transport, thus reducing the 
effectiveness of the heat pipe. Since the heat transfer 
method in the rotating heat pipe depends strongly on the 
latent energy of the working fluid, the rate of heat trans- 


fer in the heat pipe can be expressed as: 
О. = m h (I-1) 


where m, = vapor mass flow rate in lbm/hr. From the con- 


tinuity equation, 


т = р 0. А (1-2) 





where 


U ә 521507 о: the Vapor in ft/sec, and 


A = cross sectional area for the vapor flow in СО. 


The heat transfer rate becomes 


о. uM p A Ше les 
ша the vapor velocity is considered to be sonic, 
Ше = с = /ggKRT (I-4) 
where 
C Е = оса in ft/sec 
gg = зла attonalzconstamer= 32.1739 ft-lbm/lbf-sec^ 
k = ratio of specific heats 
R = gas constant in ft-lbf/lbm R, and 
T = absolute temperature in ?R. 


mee Boiling Limit 
It has been postulated by Kutateladze [2] that the 
transition from nucleate to film boiling is totally a 
hydrodynamic process. With this postulation, he determined 
the theoretical formula for predicting the burnout heat 


EM Ux: 


KG 





where 


о. 


Ad 


= KYO, AL he fo gloe - 0,77 (1-5) 


constant value 


density of the vapor in a с 


heat transfer area in the boiler in Е 
latent heat of vaporization in Btu/lbm 
surface tension in l1bf/fte 

acceleration of gravity in ft/hr* 


density of fluid in Ше а 


The experimental data obtained by Kutateladze suggested a 


ше Ісог К іп the range of 0.13 to 0.19. 


3. Entrainment Limit 


The flooding constraint in a wickless heat pipe was 


examined by Sakhuja [3], who assumed that: 


a) 


b) 


c) 


d) 


The vapor and liquid are exposed to each other 
in a counterflow movement. 

The counterflowing vapor tends to retard the 
ада morf liquid through shear stress. 
The liquid film remains stable and smooth, and 
the shear stress is usually small. 

The buoyancy forces resulting from the density 
difference between the vapor and liquid are the 


means for maintaining the counterflow. 


Ши 





e) The buoyancy forces are balanced by dissipative 
effects which are proportional to the momentum 
fluxes of vapor and liquid streams. 


ES resulting correlation for flooding is 


2 = 
- 


Е: fg 2 
E К eee) P 
where 
Е = heat transfer rate in Btu/hr 

A, = flow area in ft^ 

C = ОЕ сова, 10. 725 Бок tube with 
Sharp edged flange 

E = latent heat of vaporization in Btu/lbm 

g = acceleration due to gravity in сь“ 

р = inside diameter of heat pipe in ft 

оғ = density of the fluid in ony ee 

E = density of the vapor in пе деи 


ШО Condensing Е 


Ballback [1] determined the condensation solution 
for a rotating heat pipe by modeling the condenser section 
of a rotating heat pipe as a rotating truncated cone, from 


which the following condensation limit was developed: 


a 


О A re -т 1? 2 = 52 
E = 145 Ede Teg s w) «ЇВ. Е Lsin Ф] о - В. ) 
НЕ 311 ф 
(1-7) 
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where 


о. = total heat transfer rate in 


k = thermal conductivity of the 


in Btu/hr-ft-F 


ВЕНЕ 


condensate film 


Ре Б септесу” ol ¿lud in lbm/ft 

Ш = angular velocity in 1/hr 

ЕЕЕ = latent heat of vaporization in Btu/lbm 

T, = saturation temperature in F 

T = inside wall temperature in F 

Hr ЕИ seosıcev Of fluid in T5m/ft-hr 

ф = half cone angle in degrees 

В. = minimum wall radius in ft 

L Salon ehe wall of the condenser in ft 


The condensing limit equation 1S a function of the geometry 


and speed of the rotating heat pipe, 


 Беггісс of the working fluid. 


and the physical 


Tantrakul [4] calculated these limitations on a 


heat pipe with specific physical characteristics as shown 


in Table 1, with the results shown in Figure 2. 


ТАВЬЕ 1 


Specification of a Typical Rotating Heat Pipe 


Length 

Minimum diameter 
Wall thickness 
Шеста а ШЕ апоше 


Rotating speed 


Ще, 


14.000 


22 1101929, 
0125 
1.008 


2700 


inches 
inches 
inches 
degree 


RPM 





| 2577 
Saturation Pressure, ЕД Ше ба ту 









AAA 7,2 
8x1 06 
В8хі 0 + 
6 6 
ц ц 
SONIC LIMIT 
5 2 
| 
Е | 
3 1100 6 1409 
р BOILING LIMIT 
8 
8 
c 6 1“ 
9 
4 
S 4 
4 
> 
^ ENTRAINMENT LIMIT 
© 2 
a 2 
54 
ЕН 
1) 
9 5 
ES i05 1x10 
8 8 
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4 ц 
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Saturation Temperature То! deg F 


Figure 2 Operating Limits of a Typical, 
Water-filled, Rotating Heat Pipe 
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Obviously from the results in Figure 2, the condensing limit 
line is the predominant limitation for the amount of heat 
that can be transferred from the heat pipe. However, the 
other limitations may become important as the heat pipe 
geometry and operating conditions are varied. 
lu Eorder to augmento the heat transfer capacity of 
the heat pipe, the current efforts are aimed at raising the 
condensing limit line which may be accomplished by: 
a) a high value of cone angle, to increase the 
driving force. 
b) some type of promoter of dropwise condensation 
to increase the value of the convective heat 
transfer coefficient, h. 
c) use of an internally finned condenser to increase 
the inner wall surface area and the value of 
h, since the presence of the fin will decrease 
the condensate film thickness. 
Byemachining Zins іп the inner wall of the condenser, 


a Significant improvement of heat transfer will be realized. 


ПН МА YSIS OF THE INTERNALLY FINNED ROTATING HEAT PIPE 

A disadvantage of improving the performance of a rotating 
heat pipe with fins machined in the inner wall is that an 
increase in the wall thermal resistance results due to the 
increase in wall thickness brought on by the fins. By using 
a material with a high value of thermal conductivity, this 


ес can be minimized. 


DU 





Figure 3 





Internally Finned Condenser 
Geometry Showing Fins, Troughs 
and Lines of Symmetry 
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Figure 4 Internally Finned Condenser Geometry 
Considered in Analysis of Schafer [5] 
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Schafer [5] developed an analytical model for the case 
of fins with a triangular profile as shown in Figures 3 and 
4. In the analysis he assumed one dimensional heat conduc- 
tion through the wall. Corley [6] for this same case 
developed a two dimensional heat conduction model using a 
Finite Element Method, and also assumed a parabolic tem- 
perature distribution along the fin surface. His results 
indicated a significant improvement in heat transfer per- 
formance of about 75%, above that predicted by the one 
dimensional model of Schafer [5]. However, Corley [6] 
cautiously noted a probable error of 50% existed in the 
calculation of the heat transfer at the fin-apex, and conse- 
quently mentioned that there may be a total heat transfer 
error of as high as 15%. Also, he noted that his analysis 
was not valid fro a stainless steel condenser because of 
numerical Bee Tantrakul [4] modified Corley's 
computer program by increasing the number of finite elements 
in order to minimize the apex heat transfer error. His 
Bes with this modification converged with the results 


ШЕ Corley. 


ENUHESIS OBJECTIVES 
The objectives of this thesis are: 
l. To generate a general computer program ror two 
dimensional steady state Conduction Heat Transfer. 
Ш! TO use this program to study heat transfer in a 


finned rotating heat pipe. 
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за 


To conduct a parametric study for a proposed heat 


pipe design. 
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ПИ Ре ЕВЕ SOLUTION 


в eeVIEW OF THE PREVIOUS ANALYSIS 


Schafer {5] studied the one dimensional model heat 


transfer solution and Corley [6] studied the two dimensional 


model for an internally finned rotating heat pipe. Both 


used the same assumptions and boundary conditions based 


upon the analysis of Ballback [1], which are similar to 


those used in the Nusselt analysis of film condensation on 


a flat wall. The more important of those assumptions are: 

l. steady state operation 

2. film condensation, as opposed to dropwise condensation 

3. laminar condensate flow along both the fin and the 
trough 

4. static balance of forces within the condensate 

one dimensional conduction heat transfer through 
the film thickness (no convective heat transfer in 
the condensate film) 

Eno liquid - vapor interfacial shear forces 

7. no condensate subcooling 

2620 heat flux boundary conditions on both sides of 
the wall section (symmetry conditions), as shown 
а Hagures 5 and 6 

9. saturation temperature at the fin apex 

10. zero film thickness at the fin apex 

fee negligible curvature of the wall 
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Figure 5 Differential Equation and Boundary 
Conditions Considered in Analysis 
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Finite Element Geometry Con- 
Sidered in Analysis of Corley [6] 
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ELM cte development of Schatfer's finned condenser 


НЕО transfer is presented in his thesis. His results are 


used aS a comparison in this thesis. 

However, Since the same fluid flow and heat transfer 
theory applies to the two dimensional conduction model, 
Corley [6] utilized this model and further developed it 
using the Finite Element Method. He considered an internally 
finned geometry as shown in Figure 6, and used the fin 


condensate film equation derived by Schafer [5] 


3 
iz do (2) = 5—5 | 
Og u (Rg m mne) E COS COS g 
(Lia) 
where 
б = the fin condensate film thickness in ft- 
K вивезено и сома цселит су of working fluid in 
Btu/hr-ft-F 
оғ = density of working fluid in PAS 
ue E сове ОВ оската fluid in lbm/ft-hr 
ш) E uculanp velocity in l/hr 
Во Б пп radius in ft 
Ж = distance along apex of fin in ft 
Z = distance along the fin, from the apex in ft 
ф = half cone angle in degrees 
E = latent heat of vaporization in Btu/lbm 


He assumed a parabolic temperature distribution along the 


En surface, 
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Ши) = xr 22 РЕ (II 2) 


КООШ аге the convective heat transfer coefficient along 
the fin surface. He developed a solution of the steady 
state two dimensional conduction heat transfer equation: 


ук "УТ = 0 (II-3) 


with boundary conditions as shown in Figure 5 


а) saturation temperature at the apex, 
Е 
b) convective boundary condition in region [1] 
-k p sois - T ) and in region [3] 
an 1 San! 
ol 
т бу - посте 6) 
Pn a | 
с) symmetry conditions, ae res ten | 2 | sand 


region [4]. 
He used the 8 node quadratic isoparametric element in the 
mesh shown in Figure 6, and a modification of the three 
dimensional isoparametric Finite Element Method computer 
program assembled by Lew [7]. 
Applying the boundary condition a), T = I: at z = 0 


gives ti = ль and equation (II-2) becomes: 
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2 2 
Т. (2) = gu. 512 + ae (II-4) 


at’ 


ar 


Т - Т (2) шə -г,2 - 


E 1 512. 5) 


Beemesubstituting equation (II-5) into equation (II-2), 
 іпсепстасіпс the result from z = 0 to 2, gives an 


EusSsEron for 6(z): 


3 2 
2 2 
E ee о: po 
с? i 
ор w (Ro + x gon ф) Ше од ф cos а 
Ри \ 


(II-6) 


і 
Assuming one dimensional conduction heat transfer through 
the thin condensate film, the local convective heat transfer 


Beerfficient is: 


3 2 


ғ 0+ w (Rp т UE SOS $cos a 1/4 
] 
3 22 


Z 
EE. arp Re > 





Bu — (х2) 


IM 0) 


where Ke EN M s uml conduetivitvounoBtu/hr-ftzF which 
meomcdpplicable from z = 0 to z = 2* = Zo - 6*/cos a; where 
6* is the trough condensate film thickness, as shown in 


те 5. Across the width of the trough, the convective 


Neat transfer coefficient is: 


Sel 





Ё: 


Е = S*(x) (11-8) 
шизипеав transfer rate for the entire condenser is then 
determined within the computer program algorithm by dividing 
mal length into 100 increments of length Ax (see 
Figure 4), and solving for the incremental heat transfer 
rate using an iterative procedure. The incremental rates 
are then summed to yield the total heat transfer rate. 
Figure 7 shows a flow chart description of the computer 
program used in Corley's analysis [6]. Constant values of 
rj and s, are determined at each increment using Lagrangian 
interpolation fit of the temperature at the nodes 1l, 3, 
and 6 (see Figure 6) from the previous iteration. At the 
first increment, values of ry and S} are calculated from an 
initial guess of the temperature at these nodes. The average 


nodal values of h and the average of h of the 


lavg’ Азама 
lert hand boundary region associated with node 6 (see Figure 
6), are calculated using the position dependent value A(Z) 2s, 
(equation (II-7)) in a Simpson's Rule approximation. In 


the trough, h(x) нео сштатеа from equation (11-8), 


Puga 
Ша Се 5% value calculated from the previous increment, 
and this value is used as the average convective heat trans- 
fer coefficient for nodes 9 and 12, and the right hand 
boundary of node 6. The average heat transfer coefficients 
calculated for the left- and right-hand boundary region of 


node 6 are weighted by their respective region length. These 


values are then added to yield the average convective 
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i = 100 


ізі Ср Ах 


OUTPUT 







FINITE 
ELEMENTS 
SUBROUTINES 


9; , Thode 


Figure 7 Computer Program Flowchart of Two 
Damen sonalseondwetion Analysis of 
Corley [6] 


за 





| coefficient for node оба о The remaining boundary 


бауа” 
conditions are now applied, and the incremental heat trans- 
fer rate per unit of condenser length Qi; and the nodal 
temperature values are generated from the Finite Element 
Method subroutine. If the resultant heat transfer rate 
converges to within a 0.05$ change from the value of the 
previous iteration, i.e. 
an 


е; 


< 0.0005 (II-9) 


then the increment is deemed solved and Q; is set equal to 
the final iteration, 05: However, if the convergence 
criterion is not met, newly calculated nodal temperature 
values are used to establish a new set of inner wall сопуес- 
tive heat transfer coefficients, and the Finite Element 
solution is repeated until the heat transfer rate converges 
to within 0.0005 tolerance. The incremental mass flow rate 


in the trough is now calculated using the relation 
= 2 Ах | (II-10) 


where this value is added to the mass flow from the previous 
increments. A modification of the trough mass flow rate 
equation from Schafer's thesis [5] is used to calculate the 
subsequent interval's trough condensate thickness $* (х), 


with the polynomial root finder subroutine: 
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о = = 


. р a (Ro + xsin 9) 6* (x) 7 зіп ф Р 
и о 00° 00 Е+5* (<) tan a] 
f 


(TIST) 
Pe rementation is continued until the entire length of 


the condenser is transversed. neremental-heat transfer 


rates are summed, hence the total heat transfer rate: 


100 
Е ол) (11-12) 
я | 
where Nein 1s the number of fins along the inner wall 
Circumference. 


To start the incrementation, two initial guesses are 
made. These are the initial value of 98% which is taken 
from the analysis by Sparrow and Gregg for condensation on 
a rotating disk [8] and the initial values of the tempera- 
tures at the nodes 1, 3 and 60. 

Corley [6] noted that an error of 50% due to calculation 
of heat transfer rate at the apex may cause a total error 
of 15%. To minimize this error he suggested an increase in 
the number of finite elements. 

Tantrakul [4] continued Corley's study. He increased 
the number of finite elements between the apex of the fin 
and the trough as shown in Figures 8 and 9. To save com- 
puter execution time, he reduced the number of increments 
БО 25 instead of 100 increments as in Corley's calculation. 


However, even with the reduced number of increments, the 


Е 5 
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: > E 
Би 2 

Rout, ша 

uu Om 25-11. 

b = 0.025 in. 
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Figure 8 Internally Finned Condenser Finite 
Element Geometry (3 Finite Elements) 
Considered in Analysis of Tantrakul [4] 
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Figure 9 Internally Finned Condenser Finite 
Element Geometry (4 Finite Elements) 
Considered in Analysis of Tantrakul [4] 
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results of Tantrakul [4] converged favorably with Corley's 


results [6]. 


B. 


FORMULATION 


In this thesis the author used the same physical rela- 


tion and technique to calculate the heat transfer rate from 


the condenser. These are listed below: 


a) 


b) 


С) 


д) 


assume parabolic temperature distribution along the 
на boundary 

utilize equation (II-6) to determine the condensate 
thickness, 6 (2) 

шазиссшатспоп и 119), (1Г-10), апа (ТІ-11) to 
calculate the mass flow rate in the trough, and 
determine the subsequent condensate thickness. 
iterate to arrive at the desired solution with the 
convergence criterion 0.01, since the results are 


converged at this value (see Table 2). 


For the Finite Element Method solution, the author used a 


linear triangular finite element model as shown in Figures 


Ши пкоцай 12. 


the 


ШОО зу [6] solved the field equation for the fin with 


following boundary conditions: 


a) 


b) 


E) 


temperature at the apex is equal to the saturation 
прека ке 

convective boundary conditions in the inside and 
outside of the condenser 


symmetry condition on the sides. 
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ШИЕ first boundary condition indicates no heat transfer 
(Пе apex and a value of infinity for the convective heat 
transfer coefficient h. Since the value of h in reality 
never equals infinity, the author let the value of tempera- 
ture at the apex float, which is in agreement with the 
theoretical analysis for the extended surface in a one 
dimensional case. Only for a fin with a small angle or a 
very long length can the temperature at the apex become 
equal to the environment temperature (saturation temperature 
in this case). 

Therefore the statement of the problem for the formula- 


tion of the Finite Element Method is shown in Figure 13: 





2 2 

3 EDS A E 
EO ULIS 

9х ду 


with the boundary conditions: 


| T MN = 

Е ЙІ ІП тесісп l, -k pL hj (T р (II-14) 
Е от 

b) ІІ redgriorn Зи -k 3n = п, (E = ша) (ТТ-15) 
| н SHE 

ee in region 2 and 4, а= 0 In 


EN ding the domain into linear triangular finite elements, 


the differential equation within each element is: 


DG) NOE C, 
—— = e та) 


Эх“ ду 
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Е Di terential Equation and Boundary 
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Е: che boundary conditions: 





ve) 
-к'є) 3T ET uer. (11-18) 


m 


for elements with sides along the convective boundary, and 


эт (@) 


5), әп 





= 0 Cr IET) 


for others. ‘In the above equations, 


т(®) = temperature within the element 
ще) = thermal conductivity of the element; assumed 
to be constant 
n?! = convective heat transfer coefficient along 
the boundary of the element, assumed to 
be constant 
m - surrounding temperature at the boundary of 


the element. 


Define the approximate value of the temperature distribution 


Enn each element as: 


E 
(e)n Шер e е 
T у =) но Су) T. 
= (II-20) 

where 

М. = Two dimensional linear shape function, and 

Uc - nodal temperature. 
Using equation (II-20), equation (II-17) can be written in 


approximate form as 
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3 2 © 2 
(а) i А 

T. ен SE 5 (реа) 

: ax de A 9У 

121 ı=] 
“where 
R = residual. 
Applying the Galerkin criterion: 

[в м, да) = 9 (II-22) 


a (8) 


ema substituting equation (II-21) into equation (II-22), 


gives: 


2 2 
9 М, 9 N 
уут, (8) E N. Зд (ё) + т 
( BE oX (e) ду 
= А 


=ч, ад (2) 





А 


Applying Gauss' theorem to the factors on the right side 


yields the symmetric equation: 


ON. М. ЭМ. 
ET a ка + sac 
a (£) 
зт. 19) 
+ сс ЕК (11-23) 
"i 
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p 
|| 


area domain, and 


= 
|| 


length domain. 


Substituting equation (II-18) and (II-19) into the second 


БЕБИ Of Equation (II-23), gives: 


, 19) 1 (е) 
(11-24) 
where 
(e) 
MG - Т (е), for an element with its 
k з Side along the convec- 
С) tive boundary 
ф (Т у= 
ко others 


Therefore equation (II-23) can be written in the form: 


3 
эм. эм. ЭМ. ӘМ. 
(е) + пр аа (= 
a. хх + yay) GA 
i=l д (Ө) 
+ fa, &(T 9) aL?) = o (II-25) 
L (8) 
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Ene solution of 








о oN ON, б С 
о u = - 
= enon СЫ (різ) (11-26) 
де) 
where 
Е oy es У 
ЕНІ = E (корс у) 
17 ад (2) an (@) 
Б; = E rov (II-28) 
с. a x (II-29) 
ai 
J | = го, З 
k 
as shown in Figure 14. Recalling equation (II-24), for 
elements with an edge along a convective boundary, 
КО == _ „х= n (9! (мом. ас 9) 
3 і y Le) 2 
1 (е) L (8) 
це) 
а re) 
- TAE HE 
De сет LE | ) 
1 (е) 


Applying equation (II-30) for elements with the nodal point 


l and 2 located on the convective boundary, gives 
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Figure 14 Linear Triangular Finite Element Geometry 
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[N bam ат) = p (e) p te) f P1 ща 

3 т: x (e P2 

, 12) z (e) 0 
(е) m 

E (e) h үе 
ЩЕ к 9) | 02 aL 
(е) T 
(II-31) 
where 
01/05 - One dimensional linear shape function. 
EL sore, equation (II-23) can be written: 
О 
y] 
(e) 

(e) (e) п Фо (е) 
a) + те) f | | тя 
Ш L (8) 0 

| 
(е) 
| (e) О (е) 
(е) 0 
оба 
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where > A 


/ i / / 
p (e) (e) ) (e) 
Es c "AS рома e 217 423515: 
ins 228) aa) вк (©) PAS 
A / (е) (е) И: (е) (е / 
) 
ы! ре в” By 4e. n s bb. cC 
дА (2) gx 2) an (e) зк (©) mo 
А 
/ 2/2 
Ei... b,b,+c,c. b. 4c. A 
ад (е) "mo "mo | 
and 
(е) - (е) (е) | 1 | 33 
h T L 
(p; 9) БГК 2l — | n 


The element matrix equation is then assembled into a system 
matrix equation with the sequencing based upon element and 
System nodal point connectivity. 


Шен лез: transfer rate within each element was calculated 


by 
9 (9? ве) po er. u _ q 
/ о?) 
where 
D = temperature at node l 
i 5) = temperature at node 2 


oq 





ЩЕ = surrounding temperature 
nie) = convective heat transfer coefficient 
па (е) = length between nodal points 1 and 2 


ЕБ shown in Figure 15. 


ШЕ ТИЕ COMPUTER PROGRAM 
The program consists of a main program and three 
subroutines; 
a) the routine "COORD" used to define positions of 
system coordinate points 
b) the routine "FORMAF" used to formulate the Finite 
Element Method equations 
ШІ ene routine "BANDEC" ats an equation solver for a 
symmetric matrix which has been transformed into 
banded form. - 
The input data is entered with the units of inches for the 
“length, degrees for the angle, ЕН for convective 
neat transfer coefficients, and in F for the temperatures. 


Input data is arranged in 6 classifications, and is 


shown in Table 3 











(с) (е) 
GAF te) 
zur 


(e) 


Figure 15- Heat Transfer Rate of an Element, О 
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TABLE 3 


elassırilcatıion Of Input Data 


Elassification Type of input data 
I system and element connectivity 
HET condenser geometry 
mit statement of the problem 
IV convergence criterion 
V specification to determine the system 


coordinate point 
VI initial guess of the temperature at the 


nodal points along the fin and the trough 


A detailed description of each of these classifications 
follows: 

Classification I 

READ. (5, 100) NEL, NSNP, NBAN 


НОО FORMAT (315) 


Explanation: NEL = number of elements 
NSNP = number of svstem nodal points 
NBAN = system band width 


ЕБ 200) (TEL, (ICOR(IEL, 1) ,1=1,3) TEL=1,NEL) 


200 FORMAT (415) 


| 


Explanation: IEL the element number 


ВЫ, Г) system nodal point (NP) corresponding 


to NP I of element IEL 
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шосе: = numbering the element, first for all elements with 
a convective boundary, starting from the top to 
the bottom, then others 
- numbering the nodal point of an element is 
counterclockwise 
- the elements with convective. boundary have nodal 
point 1 апа 2 located on the convective boundary 
emassification II 
READ (5, 300) CL, CANGL, RBASE, R2, THICK, BFIN 
300 FORMAT ( 6G10.5) 


condenser length in inches 





Explanation: CL = 
CANGL = half cone angle in degrees 
RBASE = minimum radius of the condenser in 
inches 
N 
82 = RBASE # = 
THICK = condenser wall thickness in inches 
ВЕ - the eight of tbe fin in inches 


ПИО, 400) NDIV, NEST, NEFB, NBOTI, NBOTF 


400 FORMAT ( 515) 


Explanation: NDIV = number of increments along the length 
of the condenser 
NEST = number of the element on the right 
end of the trough 
NEFB = the element number with convective 


boundary located at the base of the 


fin 
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NBOTI = the element number with convective 
boundary located at the right hand 
of the bottom side 

NBOTF = the element number with convective 


boundary located at the left hand 
of the bottom side 
@lassiftication III 
БЕРІП (5, 500). (MRPM(T) ‚I=1,3) 
ВО БОБМАТ ( 315 ) 
Explanation: MRPM = number of RPM (w) 
ED (5, 600) (TSS(I),I=1,3) 
600 FORMAT ( 3G10.5 ) 
Explanation: TSS = saturation temperature of the 
копоти: (T aps) 
БЕРІП (5, 700) TINF 
BOOZEORMAT (G10.5) 
Explanation: IN Mo mes tae temperature  (T_) 
BED (5, 800) (HINF(I),I-21,3) 
ШО? ҒКОБМАТ ( 3С10.5 ) 
Explanation HINF = outside convective heat transfer 
coefficient (hut) 
Classification IV 
pep (5, 900) CRIT 
900 FORMAT ( G10.5) 


Explanation: шш = Соллесаепсе Criterion (0.01) 
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Slassification V 
ЕДО (5, 1000) (FANGL(I), (=1,3) 
NOOO FORMAT ( 3610.5 ) 
mexplanation: FANGL = fin half angle (a) 
ЕЕЕ (5, 1100) ТЕЕ 
[100 FORMAT ( I5 ) |, 
| Explanation: СІРЕ = n=, where nn is the number of rows 
of the upper triangular fin section 
m5, 1200) (KFIN(I),KFF(I), (=1,ТЕЕ) 
ШО FORMAT ( 1615 ) 
Explanation KFIN - the number of system nodal points 
located on the symmetric boundary 
e зпмоввозноціаю tin section, but does 
not include the system nodal points 
located at the base of the fin and 
the apex 
KFF = the number of system nodal points 
located along the fin convective 
boundary, but does not include the 
system nodal points located at 
the base of the fin and the apex 
m5, 1300) DOBF, DOTH, JTC, JLC, JINT, KT 
ESO0 FORMAT ( 2G10.5, 415) 


Explanation: DOBF - number of column within the fin 


DOTH = number of column within the trough 


3 





JTC = 


JLC = 


JINT = 


Note : aS an example for 
Figure 16 
Classification VI 
БЕСІ (5, 1700) T(IE) 
О (5, 1700) T(IG) 


ВОО FORMAT (G10.5) 


Explanation: T = 


number of the system nodal point 
located at the junction of the 
symmetry boundary and the line of 
intersection between the fin and 
the condenser wall 

the number of the system nodal 
point located at the center of 
„ы. coordinate 

the numerical difference between the 
two adjacent system nodal points 
vertically at the condenser wall 
section 

the number of rows within the wall 


section 


25 finite elements is shown in 


initial values of the temperature 
estimate at the nodal points l and 
2. for elements with convective 
boundary along the fin and the 


trough 
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ВВ TEST CASES 


Petore the Finite Element Method program (routine 


"FORMAF") was applied to the condenser problem, a series of 


test cases were analysed. These cases allowed a comparison 


to the results obtained from previous analyses and reference 


9. These cases are listed below: 


ШЕ; 


ессапасдлас block іп a uniformly convective 
environment on three sides, with a specified 
temperature on the non-convective boundary. These 
results are compared with reference 9 as shown in 
Figure 17. 

A rectangular block in a uniformly convective 
environment on the top and bottom side, insulated 

on the other sides. A comparison of these results 
with the exact values and Corley's results [6] is 
shown in Figure 18. 

ENrecrangular block with triangular fin in a uniformly 
convective environment along the fin and bottom side, 
insulated on the wall. A comparison of the results 


with Corley's results is shown in Figure 19. 


The results show good agreement with the results of reference 


9 (case 1), and are comparable with the exact values (case 2). 


Therefore, "FORMAF" can be used for further applications. 
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h = 120 Btu/hr-ft*-°F 
T = 1500°F 
1000.0 1420.2 МЕ 150.0 
4 2 + 
(1000.0) ПОРАМИ а si (Gis oak. OF 
k = 15 Btu/hr- re - 
| h 
1000.0 1268.5 1402.5 1469.5 1.0 ft T 
(1000.0) | Gece o la) 
Wall 
= 1000°F 
| | 
1000.0 1420.2 jm 5 1501.0 
(1000.0) | (1420.2) ЛЕ РЫ)” 
NU LE EE ue 
h, T 
29894. Btu/hr 
о =! 
(29894. Btu/hr) 
Note: numbers in parenthesis refer to the results of 


reference 9 


Figure 17 
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Statement and Results of Test Case l 








ОО E 
В = 100. ее Е 
О 
Actual T = 24.1 
ESO) (23.5) (26.0) PIES (26.0) 
o З 
245.5 245 

Y. Ed | = 15. x 


| (16.2) ec) LE ZEN m 2) (16: ый ОБЕ: 


16.0 16.0 16.0 
ye (8.6) (6.4) 
p. 6 Deb 1.6 / e 


Actual T = 7. 2 F 
ODE 





ше ЖИЕ 


LOO ORO Btu/hr-ft^-deg F 


EJ 
|| 


, 7586. Btu/hħhr actual 
7552. Btu/hr 
Про ВЕЩИ Не 


И 


‘Note: Numbers in parentheses refer to results of 
Corley [6]. 


Figure 18. Statement and Results of Test Case 2 
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N hy ГОО ЕН hr- fe deg E 
COPPER | 
- 0/55 
j TE 
T = 100 
2 
h, = 1000. Btu/hr-ft^-deg F 






HEAT TRANSFER RATE (Btu/hr) 





ONE TWO TWO 
DIMENSIONAL DIMENSIONAL DIMENSIONAL 


6526. ои 6826. 
6561. 6754. б 750. 


*) Results of Corley [6] 





Figure 19 Statement and Results of Test Case 3 
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ІІІ. NUMERICAL RESULTS AND DISCUSSION 


The analysis concerns itself with the geometry of both 
 Соррег апа stainless steel condenser as tabulated in 


Table 4, using water as the working fluid. 


TABLE 4 


Exeocrfrcation of a Typical Internally 
Finned Rotating Heat Pipe 


Condenser length = o O00 inches 
Minimum radius = 057 75 inches 
Wall thickness = 0.03125 Inches 
Height of fin = 0502 5 inches 
[ Cone half angle = 1.0 degree 


all 


- 


Results are given in Tables 5 to 18 and Figures 20 to 28. 

It can be seen from the data that the solution of the Finite 

Element Method converges. This can be seen by examining 

data of temperature distribution within the fins of 25, 

40, and 108 linear triangular finite elements (Tables 7 to 

9), for the same location (nodes 10, 15, and 28 respectively). 
The heat transfer rate solution is generally affected 

in the same manner as the analysis of Schafer [5] when the 


outside convective heat transfer coefficient (h J)e rotational 


QuE 


) 


speed (w), fin angle (2a), and thermal conductivity (k 44, 
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TABLE 9 


Theoretical Heat Transfer Rate on a Stainless Steel Finned 
Condenser at 1000 RPM, for nout = LOCO, Btu/nr- fto deg.-F 


half of се 
angle fins POM 108 EL 






















































10 276 | 8948.2 8944.1 8893.9 
23764.0 23727.0 23607.0 
38557.0 38492.0 38286.0 

30 84 8515.5 8494.4 8461.8 
22511.0 - 22445.0 22354.0 
36467.0 36355.0 36204 .0 

50 40 8186.3 8164.6 8141.6 









2u5 740 
34916.0 


2750820 
34805.0 


21020 
34693.0 
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TABLE 10 


M@meoretical Heat Transfer Rate on a Stainless Steel Finned 





Condenser at 3000 RPM, for n = 1000. Btu/hr-tt*-deg.F 
P A Dr QUSEEU M) 

O 
angle fins => Ова БЕ 40 EL 














































10 276 100. | 9096.4 9100.4 9058.1 
150. | 24193.0 24200.0 24098.0 
200. | 39277.0 39284.0 39123.0 

30 - 84 100. | 8818.4 ° 8810.3 8786.0 
150. |23389.0 ОО 23301.0 
200. | 37934.0 37900.0 37781.0 

50. 40 100. | 8584.1 8580.7 8561.1 
ИН! 22690.0 22642.0 
200. | 36823.0 36768.0 36689 .0 
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TABLE 11 


Theoretical Heat Transfer Rate on a Stainless Steel Finned 


Condenser at 1000 RPM, for h__, = 5000. Btu/hr-ft*-deg.F 


sat Q(Btu/hr) 
half of (OF) 
angle fins Zo EL | 40 EL 













10 22 100. _ 26899 .0 26775.0 26308.0 
150. 70520.0 68673.0 
200. 113950.0 113300.0 110780.0 

30 84 100. да ДО ЩО 22540.0 22286.0 
150. 58879.0 58288.0 57582 ..0 
200.. 94734.0 93748.0 92581.0 

50 40 100. 20426.0 20253.0 20086 .0 
150. 52480.0. 51978 .0 51513.0 
200. 84275.0 83443.0 82676.0 
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TABLE 12 


Theoretical Heat Transfer Rate on a Stainless Steel Finned 
Gondenser at 3000 RPM, for hg. 7 5000. Btu/hr-ft?-deg.F 








































































"in 
заїї 
angle 
10 276 О 28.1010: 0 28263.0 
O 75734.0 74492.0 
722610920 1225200 72057020 
30 34 25604.0 25334 .0 25151021: 8, 
66822.0 66734.0 65951.0 
107830.0 107800.0 106260.0 
50 40 231673 50 23550 .0 2 9h E959. 








60726.0 
97803.0 


611270 
98464.0 


61488.0 
99069.0 
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TABLE 13 


pemperature Distribution Within the Fin at the 
1051 Increment. 
(Copper Condenser with 25 Finite Elements 


deg 
1000 Btu/hr-ft^-deg F RPM = 


= 496 97.2 Beu/nr 





= 100°F Q 


T3 





TABLE 14 


Temperature Distribution Within the Fin at the 
10th Increment 
(Copper Condenser with 40 Finite Elements) 


L 
2 
3 
4 
5 
6 
7 
8 
9 
1 


deg О 
= 1000 зу Два чес Е ВРМ = 1000 
100965 9687.2 Btu/hr 
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TABLE 15 


Temperature Distribution Within the Fin at the 
10th Increment 
(Copper Condenser with 108 Finite Elements 





99.04 24 98.39 50 97 90 
98.91 у 98.05 
97.96 
97235 
97.84 
97-.97 


92.36 
91.95 


рн о oa 4^ со NY Hn 










а deg о = 

hout = 1000 Btu/hr=£tó-deg Е RPM = 1000 

T = 100°Е О = 9678.9 Btu/hr 
sat 
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TABLE 16 


emperature Distribution Within the Fin at the 
10th Increment 
(Stainless Steel Condenser with 25 Finite Elements) 


Dep [= [a [= [oon 
99.79 15 93,7 


99.08 





юк о со ы соу Ui > (9 + 








50 deg о 
= 1000 But/hr-ft^-deg Е RPM = 1000 
200. 7 Био, ВИДЕ 
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TABLE 17 


Temperature Distribution Within the Fin at the 
10th Increment 
(Stainless Steel Condenser with 40 Finite Elements) 


Се ен ЕЕ НИ 
99.82 15 96.70 29 


о ; 30 





к 0 ы су (л ы U N ғ 






Е deg 7 _ = 70°F 
houg = 1000 Btu/hr-ft*-deg F ^ RPM = 1000 
T | - 100?F DENN ERE 





sat 
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TABLE 18 


Bemperature Distribution Within the Fin at the 
10th Increment 
Stainless Steel Condenser with 108 Finite Elements) 


ете е е ре ра 
| 99.84 24 Don 50 94.46 
29.50 26 3 94,33 
94.64 
94.29 
В о 
33,16 


96287 
22:84 






HM oo N боб wv мо ыны 





deg T = 70°F 
1000 Btu/hr-£t--deg Е RPM = 1000 
100543 3154 | -бавЕш/Аг 
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а = 30° 
= соррег 
- - = - = stainless steel 
1 
160 
b. 
E 129 
5 
42 
= 
er 
г | 
е = 
© 80 a 
40 
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Saturation Temperature (°F) 


5000 ИН Ее — deg Е 


x Aout в 
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. = = ze 
2 E 1500 Btu/hr-ft есе 
I = 1000 Btu/hr-ft^- deg Es 
QUE 


Figure 20 Heat Transfer Rate (0) ої Internally 
Finned Condenser at a posts cu vValue 
fen vs. Saturation Temperature of 


/ 
Sue Working Fluid 
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RPM = 1000 2 
ШЕ е = 1000 Btu/hr-ft -deg F 
=- = copper 
- - - = stainless steel 
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Figure 
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302 
50° 


za 


Е c O = a Pe а і 
150 200 
Saturation Temperature (°F) 


Heat Transfer Rate (Q) of Internally Finned 
Condenser at a Particular Value of Fin Half 
Angle vs. Saturation Temperature of Working Е 19359 
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Dese ZEHDUOOSOBtu/hr-rt =а едЕ 4. uc 13000. Btu/hr-ft -—degF 
bs t Е: | _£+2_ 

De = 5000. Btu/hr-ft -degF 5. E im0COne seu, hr-Lt =degr 
n 2 

Desi = 9000. Btu/hr-ft -degF 


Figure 22 Heat Transfer Rate (Q) of Internally 
Finned Copper Condenser vs. RPM 
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Figure 23 Heat Transfer Rate (Q) of Internally Finned 
Stainless Steel Condenser vs. RPM 
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Figure 24 Thermal Resistance of Internally Finned 


Copper Condenser vs. Е 
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Figure 25 Thermal Resistance of Internally Finned 


Stainless Steel Condenser vs. nest 
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COPPER CONDENSER 








RPM = 1000 
| 0 
cat = 100.0 ÔF 
70.0 °F 





?- -г-т- 
о 5 10 
b 
a 
Іі 540.01 ы 548 
а а 
ES o1 ША 16-4. 15 
а а 
3. oe 0-5 6. E z 2 
a a 


4” 





heat transfer rate of 
internally finned con- 


denser 


: heat transfer rate of 


smooth condenser 


Figure 27 Comparison of Heat Transfer Rate (07/05) 
VS. 2 ato = r000: 
a © 


ULC 
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Btu/hr-ft^-degF 





ту 
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4 
5 
m 
3 
COPPER -CCNDENSER 
RPM 2090 
e T = 100 °F 
sat 
To = 70 OF 
БЕРГ -т——т— 71 
0 5 19 b 18 7 20 
а 
1. ENE 0.6] Е Q : heat t f t 
= : Е : г! hea rans er rate gf 
5 (с o internally finned con- 
=a 0,1 — яв ],5 
а а denser, 
Bee с c 
= 0.5 5 2% Re : heat transfer rate of 
smooth condenser, 
Figure 28 Comparison of Heat Transfer Rate (9-/0-) 
vs. 2 ath = 5000. Btu/hr-ft^-degF 
a HULE 
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were varied. In Figure 20, the influence of the outside 


convective heat transfer coefficient, h 1s shown both 


Quer 
for copper and stainless steel. Figure 21 illustrates how 
the rate of heat transfer was affected when the fin half 
Нисе саз varied. In addition, the heat transfer rate was 
found not to vary significantly with rotational speed when 
the heat pipe operated with low values of outside convective 


heat transfer coefficient, h However, for high values 


out 


of h the effect of rotational speed cannot be ignored 


sut ' 


— 


(see Figures 22 and 23). 
The small variance of the heat transfer rate with rota- 
tional speed is due to the dominancy of the outside thermal 


resistance, R such that changes in the condensate film 


sur! 
thickness, 6, will not Significantly affect the total thermal 
resistance as shown in Figures 24 and 25. The heat transfer 
rate will have a large variance when the heat pipe is operated 
with a high value of hout' 

From the results of the parametric study, as shown in 
Figures 27 and 28, the heat transfer rate continuously in- 
creases as the fin half angle decreases (b/a increases). 
Since the total number of the fins in the condenser is 
larger with a small half angle fin, this results in an in- 
crease in the heat transfer rate. However, the increase 
Nmoniy slight when the fin half angle less than 11 degrees 
(b/a greater than 5). 


The study shows the enhancement in heat transfer is 


larger when operating with a high value of Aout or when 


88 


m 





E. is small. Therefore the design operating with a 


Dach value of AS using both internal and external fins 


ut" 


is recommended. 
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ПОСНА О ТОМО 


The Finite Element Method works and converges. 
The designer has a choice of: 

- material 

- number of fins and fin angle 

E] MA 

= Aout 

The heat transfer rate increases continuously 
with decreasing fin angle; however, for half 


angle less than ll degrees, the increase is 


only slight. 
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V. RECOMMENDATIONS 


1. Manufacture internally finned heat pipe and test. 
2. Manufacture internally and externally finned heat 


pipe and test. 


S 
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APPENDIX 


SENEUTER PROGRAM LISTING 


OE TIMENSICHUSL HEAT CONDUCTION FINITE ELEMENT ANALYSIS 
Meee leniNG HEAT PEPE . USING TPIANGULAR ELEMENT MODEL 
ПАШЕ Py MAJOR IGHATIUS.S.PURNOMO IN JUNE , 1973 


X 3 06 Эс 5 5€ 3€ X 


ШЕГІСІТ РЕНІГ(2(а-Н,0-2) 
ПЕТИ 21200) EPS(2009),H2(200).XC0F (5 ),COF (5). ROOTF(4),RUOTÍC(A 


IED, GE 200).0MDOTL200),UF (2900), CF (200),RHOF (262), DEL(200),CWaA 
x09).AMTOT(209),R(2080 ).01N2 (200), T3(200),TT(209)3,TB 0 (299), TE (299 )5 


ee Geo). NCCCOO0) 

ПА М! Т55Е3 3, РАМОС. (3), НТМЕ (3), ИКРМСЗ) 
EENUUSBZABPOL^ADOBF,DOTH,KFINGIOO ,KFF(CSOD,IFF.JTC, JLC, JINT,KT 

В ВОМУМАЕО/ Е (209, 50),.2Е (290,1 ),Н‹е99), 151209). ТЗАТ ‚ СК, МЕ, ММ? . Нл 


XN. ICOR(220.3) 
EON PCRD/ Xí 209), Y (200), EZEPO, BUIN, THICK. TAaLFA. APS 
BENSENCAT,B1,ZZ220*(-1.0DOfCA1XZZXtE5/3.0DO*B1xXZZtfX27/2a 00082) 


aE RENT CONNECTIVITIES 
ADS, 1890) NEL., HSHP, MBAN 
ESTIS? 
E UECS,150) NEL, NSKE. NTAN 
meri 7ex. “SO.OF .ELEMENTS=” , IS, 104. NO.OF SYSTEM H.P.*^.IS,10X.^ 


XO OF BANDED=",15) 


о 
INO 





2 


м" 


39 


РЕСЕ) KIEL,CICORCIEL,T2, Ie1,3). IEL-1.NEL) 
FOERFIATLAIS) 

LETTEL. ZSO) 

EE SS ECENENT ,16Х."НРІ",14Х,'НРО',15Х, "НР3') 
ESTO TEL,  ICOR(IEL,.I).]1=1.39, IEL=1,. NEL) 
шасито, 12х. 15, 12<, 15. 1е<,15) 


ШЕ СОНОЕМЗЕК СЕСМЕТКУ 

Сс THE CONLENSEPR СЕНОТ 

aae 15 THE CORE HALF AIGLE 

RBASE IS THE CONDENSER RADIUS AT THE EASE 

ит THE WALL THICKHESS OF THE CONDENSER 


READIS.200) CL,CANGL,RBASE,R2. THICK. BFIN 
FORMAT(6G10 5) 

ШЕСЕ 12 ога 

R2=R2-12 CDA 

RPASE=RPBASE/12 ADO 

BUIN=EFIN/12 508 


ЖВІШІ) TS THE NUMBER OF INCREMENT 
ШЕСТУ Wide ELEMENT NUMBER AT THE END OF THE TROUGH 
hems 1S THE ELEMENT HUMBER AT THE BASE GF THE FIN 
ПИ То Tre FIRST ELENENT AT THE BOTTOM SIDE 
МЕЗТЕ IS THE LAST ELEMENT AT THE BOTTOM SIDE 


READ(S. 499) MDIV.NEST,NEFB.NBOTI,NBOTF 
МАТ (515) 

DIV=DFLOATCNDIV > 

Be 141592652509 73900 

PHI =2 OD@XCANGLAPI/’SES 909 

eee =DSINC PHI ) 

EXUI"DCOSCPHI) 


^. TPHI*DTAN(CPHI) 


эз 





DELX=CL-DÍIU 

CEASE*+ 2 ОГӘЖРІЗЕКБазЕ 

БЕҮХІТАРЕНСЕЗСІ ЖТРНІ 

ШЕХІТ-О ОротРІЖРЕХІТ 

THICK STHIVH 1: . 909 

ПИ о, 7oU CL, URHGL, Re. THICK .REASE, BUIN 

50 FORMAT? 9.7759. ’CONDENSERA~S LENGTH’, Ete .S.30X%, “HALF CONE ANGLE~ 


EN 577-54, 'R2-",Eig.5,45^4,'THICK*",E12.5,7,5X., 'R-BaSE-',E1g.S5,^ 
ПЕ! = ^,Е13 5) 
DATA FOR RUNNING 


| РЕар(5.500) (MRPMCTID.I-1,3) 
ша ЕОКМАТ: 315) - 
| READ(S.600) (TSS(T),1#1.3) 
ED ЕОЕМат с 3019 5) 
ВЕС. 7099) ТН: 
39 FORMATIGIO.S) 
NEADCS, Z200) (HINFCIO.I*1,23) 
20 FORMAT(3G19.5) 


fee CONVERGENCE CRITERIAN 


READ(S.900) CRIT 
00 FORMAT(G19 9) 


INTERNALLY FIN GEOMETRY 


$ НЕВОО5,10060) C(FANGL(I).I=-1,3) 
PO FORMAT (310.5) 

НЕНОС5.1109) IFF 
20 FORMATCIS) 


94 





210 


32 
| ННЕ=НЕРВ/ 2 


ша оо ест) (СКРТИАСТО,КРРСТО, 1.1 .ТЕЕ) 
цшсмат (1615 ) 

Оо, 109 3р7ВЕ , ООТН. УТС, Л.С, SINT. KT 
А сес[0.5,415) 


ВИ = НЕОТТ+ЕМНБОТЕ-НЕОТТ }/Р 


Ж 


ЕЕ ЗНЕОТЕ + 1 
ПО 15350) ICORCNBOTI.€), ICORCNEFB,.1).ICOR(NTHM.2).ICORCNEST. 1^ 


ECORIMBOTF .1) 


BE IU EZZ5X,'TIB-'/.15,10X, 'TT*", 15, /. EX, 'TEM*s 7, IS, 10X, ' TES", IS. ^ 


| 


x 


me TES^.IS) 

ПО 14020 КІн«1і,3 
ALFASFFANGL(KIA же 0DOXPI/3€690 9109 
SALFA=DSINERLFA) 
CALFA=DCOS[(ALFA ) 
TALFA2DTANLALFA ) 

EZERO=2 GOCABUINEATALFA 
EPSO=EZERO й 
MFIN=CEASE-(EZERO+EPSO ) 
EPSEX=(CEXIT-(NFIMNXEZERO))/NFÍN 
NAS TEPSEX-EPSO3)/DIV 
ZZERO=BUIN/CALFA 


. ZA=9.0D0 


| BOUNDARY CONDITIONS AND TEMPERATURE *S ESTIMATE 
ALONG THE FIN BOUNDARY 


M" 


ро 1459 NTINT=MBOTI.N20TF 


Ша 1 SCNTINEG D=TING 


DO 1500 NNT«NSF,NEL 
TS(NNT)s0 2509 \ 


9 HONNTI=G 9D9 
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DO 1600 IGT=1.NE=T 
IE- ICJRLIGT,2) 

EENEUS 1700) T(IED 

IGs ICOR(NEST, 1) 

ESSO 1700) ТСТС) 
FORMAT(G10 5) 

ГО 1490 KIB*1,3 
ПЕРМ = МРЕМСКІР ) 
OMEGAsNRPMYr2 VUDOXPIXGO BDO 
DO 1400 KIC=1,3 

DO 1200 KL=HBOTI,NBOTF 
EE) 2 HTTIF(KIC) 


HIFN=HINFCKIC) 


Boa КІр=1,3 

ВОНО TSS(KID> : 

DO 1300 HSRTs1,NEST 
TSINSAT 22 TSAT 
POBID>=(TSAT+TÍINF 3)/2.0D0 
DEL(1)=3 YOH0O6TS52D9 
QT-9 ODO 

QTOTsO ODoO 

DMTOT*O ODS 

NK =NDIVU+1 

De72069 NI=1.NK 
RCNT)»R2+NIXDELXXSPHI 
В: 2=ЕРСО+НТЖВЕТЯ 
ВЕРУ СМ] } 


NODAL POINTS COORDINATE 
CALL COORD 
2(1)*2А 
DO 2109 IZEU=1,NEFB 
NA*ICORCIZEL,1) 
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то 


ow 


90 


BEXTICOP:IZEL,23) 

AXERXUCHAO-SCHP) 

ТЕ < г на )-У (МБ ) 

ЕТ7 < рзсРТ. ЕХО ТЕЖАЕВЪ! 
НЕ +] )=2. 1220 )+Е(2 

МИР ГТСОРГГІНБ. 1))-Х(ІСОРС(1,2)) 
шаси ите АНЕ, 11т)-У(СТСОВС1, с) ) 
ZbESDSCOPTCOXZBEXOQT-YZBKXQ) 

DORE ZERO 

IM=1 


PARABOLIC TEMPERHTURE DISTRIBUTION ALONG TRE FIN 
BOUNDARY USING LAGRANGE INTERPOLATION 


Seen © [COR(1,.2)) 
ОО ЗТСІСОНСІЯНЕ.1 1) 
TP3sT(CICOR(CHEFB,1)) 
AP1=TP1-(2BX20) 

eee etre (ZEA CB-2C)) 

Aap3s TP3/CZCXCZC-ZB)) 
Peat ZF +2C KAP 

ВРез-2СхаРӘ 

ВРЗ=-28ЖАРЗ 

A1L=AP1+4OP2+AP3 

B1=BP1+EP2+BP3 

TCsO сера 

DO 2200 NY=1. NEST 

SS E TCICOR (NE, 20) 

AY2DFLOAT  CNY+1 ) 

TF =(TC+TCICORONY., 1) tAYXTSONY D6 Ce. ODOKAY ) 


EBETD-FLUID PROPERTIES 


НЕС «10597 .2р9-2.69187500Х75 01) 
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RHOF (NI )=62.77T4D0-0 0025563810ЖТЕ-0. #9005357209 ТЕ жа 
CFUNI)=0.3034D0+0 000732927DOXTF-0 00000147321D9XTFXx2 
UF(MI)=0.001327D0-0 00001 44B39DOXTF+0.0008096631253DOXTFXX2-0.00%5 


жодовезто5бездпохтежжз 


UF CHI )2360Q0tUF CHI?) 

МОК о-=51. гге0ро-0.огагзарохтозор 

CKsCUCHI) 

COHNSTsSRHOF GNI Of *2xOongEGAKxaxHFGECPHIXCRLEAXRCHI) 


AVERAGE ELEMENT CONVECTIVE COEFFICIENT ALONG THRE 
FIN BOUNDARY 


ESUSBESZZERO-DEL(CNI)^Z/CALFA 
AZZ=DELCNII/SALFA 

227Z9TAR 

AZS=DABS( 4*CF (NI XUF (NI DXHDEN(CAL, BL, 2Z CONST D448. 25D0 
HAC=9 009 

DO 2303 IEL*1,NEFB 

Ew CIEL) 

BESSECIELT1) 

BEEUZSTAR LE.BZO) GO TO 2400 
GO TO 2500 

IFCHAC NE @ ODO) GO Т0 2690 
BZ=ZSTAR 

ШЕСТЕТ NE 1) GO TO 2769 
AK=(BZ-AZ)/S.0D9 

LLAK 

GO TO 2809 

AK= (B2-AZ)/4.0DO 

2282 

ZEL=AXAK 

DO 2920 NH=1,5 

HZCNH I=DABSCCEFLNIIKKIXKCONST/CAZUFCNIIFIHDENGAL,B1,ZZI I KO. 2SDO 
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20 


209 


120 


20 


20 


1139 


ЖА?) 


22-СссжғаАак 

CONH2AK CHE CL 4442 02 +CKHZ 03 )4+4KHZ 64) 4HZ65 ) 7 (3Х2Е  ) 
IF CESTAR .EQ.BZ) GO TO 3000 
НСІЕІ,) аз СОМН 

бо ТО 23409 

ALZzCSTAR 
HAZ=ENNHKIAZ-ZUIEL)) 
DELA=-AZS 

рими (ТЕ +1) 
DELB=(BZ-ZSTAR )¥AZZ/ CZZERO-ZSTAR ) 
Meee *OPelLAtLELB Ive ODO 

HAC? (EZ-AZ ECF CNI)/DELZ 
НСТЕТ) - CHAZ+HACITCRZ-ZUIEL)I) 
GO To азоо 

воза СІ) 

DELA=DELE 

HaZ=9.0DU 

GO TO 3100 

CONTINE 

NETI HEFB+1 

293209 IEL=NETI. NEST 

БЕ СЕ СГ) УБЕБ СМТ ) 


ет INTO THE FINITE ELEMENT SOLUTION 


CALL FORMAF 
CALL BANDECCNSNP.NBAN. 1) 


ШЕЕ ТЕПРЕРАТЦОКЕ DISTRIBUTION 
ШО 3300 NT=1.NSNP 
MINT) =E (NT, 1) 


ВОН )«ТСІСОНСМЧВОТІ. 20) 
femetyl JeT CICORCNEFB.1)) 
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3-30 


3:30 


3120 


330 


129 


| 
| 


meron. =) CICGRCNTSI.2)) 
TECNI ет (ІСОЮІГІЄ ОТ, 1)) 
Peet TE ICORCNBOTF,1)) 
TTS=0 еро 

Du 3400 N5=1.NSNP 
TTS=TTS+TINS) 
PN=DELOATIONS } 
TS9LID-TTS/PH 


с ат ТНЕ БОТТОМ 5105 


ОРІзо ODA 

DO 3509 IBEL=NBOTI,NEOTF 

NKA=ICORFIFEL, 1) 

HKB=ICORCIBEL, 2) 

XB=X CNA I-XCNKB } 

YEsViHER)-YCNEB ) 

ELBEDSGRT (XBAXC+YB4F Ke ) 

Gere ee i+(TCHKA)+TCHKE )D-2*TSCIBEL) )SELBXHCIBEL2/2.0D0 
ОБІАНІ )«ОБІЖОРЕЇХ 


Nu -OUrTON UNTIL CONVERGENCE CRITERIA IS МЕТ 


EESCIM-EO 15 GO TO 52609 
GJ>0B1 

GO TO 3700 

GI=GB1l 

и >= 

Ge TO 2 
аОо=разосој-01)70Ј 

МАСКЕ СКІТ) СО ТО 3809 
QI*GJ 

GO TO 2 

DMDOTCNI )#=2*QBIXDELX/HFG 
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DMNTOCT =DATOT + OMIDOT tnt ) 


Сі =КНОЕ МІ Sk KOMEGAKROXRONI D7 CSKUF CNI) ) 


ВОН N=-DMTOT 
ACOF (2120 ODO 
ALOF CS )2@ ODO 

XCOF ( 4 )C1XYEPS CHI) 
есеј С1хТАСРАа 
[12 4 


ПОНТ ГР ТСАСОР . СОР, 1], КООТК, КООТТ,ТЕР) 

дођро AMD ROSTR(1) бт е. 
09STO AND. ROOTR(2) GT 0. 
COSDO UND ROSTRO) GT.O 
BESDO.AND ROOTR(4).GT.O. 


ШЕ(РПГОТРС1) LT. 
ШЕ (РОСТРСО) ІТ 
КОО ТЕЗ) LT 
BESUROOTRIAPO.LT. 
URITE(&6.4300) 
ent 710%, “CRASH, CRASH, CRASH’ ) 
ИРТТЕСВ, 4350) (KROOTR(CIO,I*1,4) 
ҒОРМАТ(//5Х.,.4(Е120.7.3Х)) 

GO TO 6108 


колою ри» 


THE CONDENSATE THICKNESS 


DELCNI+1 )=ROOTR(1) 

GO TO 4429 

БЕГСИІЗІ )"ЕООТЕСО) 

GO TO 4400 

DELCNI+1 )#ROOTR( 3) 

СО ТО 4469 
PECONI = ROOTR( 4) 
QEL=3 009 

IF(NI NE 10) GO TO 4590 
WRITE (6. 45090 ) 


000) GN 
оре) GO 
209) GO 
орао) GO 


РР | 775X,7°NP’ 6м, х", Ех, "ү, ах, Тт” ) 


DO 4722 НР=1, М5ЧР 


URITEC6,4809) НР, Ж(МР ). (МР), Т(МР) 


TOL 


TO 
TO 
ше 


то 


3392 
4209 


AIAN 


та” 


4299 





С 


ipo 


9 


509 


(29 


Meera (7en, 13, 3. 2°F13.6,3%)) 

WRITE (5,4209 ) 

ЕСИР EL. 3x, H, 114., ' EL-LEHNGTH^,15X,/'Q-EL'*) 
DO Босе EXL=1.,MBOTF 

ИКАЗТСОР(СККТ, 1) 

NKY=ICOR(KKL. 2) 

ER A ) 

РУСО )- Y CHKY ) 

EXYsDSOQRT(!XPXx*2Q-vPtta) 

GEP DAES © CTCHKXI+TCNKY )-AKTSCRKL ) I)XEXYXHCKKL)7/2.0D2) 
GCEPsQEP XDELX 

КОО EJCO) KKL,H(CKKL :, EXY.QEP 

EE “= х . 12, 3х,Е12.5,3Х,Е18,5.10Х ,Е12 5) 

ЕТО. 5200) СРТ 

FORMAT (-2X, * CONUERGENCE CRITERIAN=**,E15.8) 


EE ROm THE TOP SIDE 


En 500 IGEL*1,NEST 

KAZ ICORCIGEL.1) 

KB=ICORC(CIGEL.2) 

ХОЕТ «ХОКА?:-ХОКВ ) 

ШШЕБЕУ(КВ )-У(КА) 

ELIT*DSGRT (CXGELXX2O *vGELX X2) 
QEL*QEL-*C(OXTSCIGELO-TCKA0-TCKBOOXELIAXHCIGELO/2.O0DO 
CONT INUE 

QINC (NI )=GELXDELX 

AMTOTCNI )=DMTOT 

GET *QELXDELXXNF INXO 

QT =QT+GET 

QAsGBIXDELAXNF INKe 

QTOT*QTCT-*GR 

CONTINUZ 

WRITE(6.S4@@) HFG.NF IN. HCNBOTF ). TSAT.NREM, QTOT, QT,FANGLCKIA) 


1518/77 





| > 
40 FORHRT(' ', 4/7, 5X, "HFG* ^, E18. 5, 7, 5X, "NUMBER OF FIMs'^,15,7, 5X, "H-O^ 


BEN EISE LZ,SX,'TSAT*'",E12.5,^/, 5X, '^RPM*^,IS,^/,5X, "^Q-BOT«*",E12.5.7^ 


ED cre =",612.5,/7,5X. “HALF-SHGLE=’,F2.3) 
WRITECG. S500) 
ЕЕ EOBlInTC'O^,ex, "^ J^, AX, 'FILM THICKHESS^.6X, 'G-IHCPEM^,6X, "haSS-TCT^ 


EN 2^4.'TT^,10X, 'TE*, EX, ТВ“) 
DO 5600 M*-1.NDIU 
| PE 555) ЧЕ. ПЕ (Ме). GBCNRD,ANTOTCNR), TIBCNR),. TTCNR),. TECNR),. ТА 


ЖЧ?) 
МИНА: °.4%,14,4X.F12.10,4%.F10.4,6%,F9 S.GX.FS.1.6X.F5.1,6%.- 
ШЕ ex,Fs 1) 

URITE(6.5200) 
bo FORMAT('0',6X. ^J'.GX, 'K-UALL' , 4X, "K-FILI^ , 3X, "DENSITY", 4X, "UISC-^ 
| 


XILII^.6X. "EPSILON^,SX, "RADIUS'/, SX, ТВМ",5Х,'О-ВОТ") 
DO S900 NG*1,NDIU.2 
00 URITE(CB, 77?) NG.CUCNGO, CF (4G), RHOF CNG) , UF CNG), EPSCNGO, RCNGO , TBI CA 


X).QINCING ) 
Р FORMAT (* °,4X.14.4X.F7.3.4X,F6.4.4%,F6.3,4%,F9.7,4X,F9.7,4%,F7 ~ 


emer > 1, 1х.,1Р1012.3) 

2Q CONTINUE 

9 STOP 

d EN D 
SUBROUTINE COORD 
ша ТОТ? REALXS(A-5H.0-2) 
COMMON/PCRD/K( 200). 7¥(2S0),EZERO, BUIN. THICK, TALFAR.APS 
memnOn7~ReOL/LORE , DCTH.KFINCSO).KFF(SO), IFF.JTC.JLC,. SINT. RT 


BOS 





514 631 


DELH=BUIN’DOEF 
AC1L)=E0 еге 

Y11)=THICK+EUIN 

=] 

ШИ 321 Ізі.ІРҒ 

ICAFKF INCI) 

ICB2KFF +I) 

CEA=DFLOATCICB-ICA) 

AN=0 009 

зе 11-Тен, ICR 
ACII)=XC1)+NAANKDELHXTALFA/CHRA 
eee = © 1 J—-NEKDELH 

ANSANT1L ADO 

Мај] 

анзо ора 

МЕРЕ ТоВ-ТІСА+1 

нене Ss = УТС, ЛЕС, ЕТ 

КО Ух» 1) 

У(Т)ес1 ODO-AN/DOTH )XTHICK 

292322 JJs1,ICD | 

СЛ tJ IXZEZERO-’ Cex (CBA . 2р0) ) 
ВТ) =) 

Шон 365 К+1, КТ 

ХО + +Коех( + 77) +КЖжаРзи се. еПожкт) 
Y(J+JJ+K)=Y(J) 

AN=AN+1 ёро 

RETURN 

END 

SUBROUTINE FORMAF 

IMPLICIT SEALX£(a-H,0-2) 

EM SHUSIOM BC35,C(32,ERnC3, 3) 
COMMON+PCRD/X(2860).,Y(289),EZERO, BUIN. THICK., TALFA. APS 
COMMON/MAFO/A(288,59).F(286,1),Ht200403.T5(2090),TSAT,CK,NEL.NSNP, Na 





BN, [CORE 280,35) 


DO 101 N=1,NSNP 
ЕСМ, 1 120.009 

DO 102 MAs1,NEAN 
ам, Ма) о. оре 

CONTINUE 

PORTOS 1221, МЕС 
ІачІСОРІТЕ.. 1) 
BESICOR:(IEL.2) 

mee lOUrRCIEL, 3) 
BC1os/".IBO-"" CIC) 
СТС -У СТА) 

BCS «"%CIRI-VCIB) 

Ree = )—х 018) 
CCaQosXCIPRO-XCIOCD 
ESSUSACIE)-XCIOR) 
EL=DSGRT(CC3)**¥2O+B13)KxX2) 
AS=DABS((B(1IXC(2)-B(2)XC(1))/2.0D0) 
HC#H( TEL)/CK 

DO 194 J21,3 

Bem CORCIEL.J) 

ШОО Об К=1.3 
KK*ICOR(CIEL,K) 
EACJ,K20s(B8CJOXBOKO-CCJOXCCOK00/7CAXAS > 
IF(HC EG 0 ОГО) GO TO 106 
HELsHCXEL/6 ODO 

ТЕСТ ЕС 3) СО ТО 106 

ТЕК ЕО 3) СО ТО 196 
ВЕС К) GO TO ic? 
EACJ,KOsEACJ,KO0-*HEL 

GO TO 106 
EQRCJ,KOoxEACJ,KOT2OXHEL 

EN ГТ JJ) GO TO 105 
NUSKK-JJj*1 
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НЕ, НУНС, НЫ )+ЕЯАС У. КЮ) 

й CONTINUE 

ШЕ СОНТІГІЛЕ 

ЖЕРЕСПСТТЕСІСІ ЗЕЕІ/2 ODO 

Peete 2 S=F CTA. LItFE 
Ри =F CIE. LI+FE 

CONT IAE, 
КЕТИР 
END 
SUBROUTINE BANDEC (NEG. MAXB,.NVUEC) 
NIE TT PEALXS(CA-H, 0-2) 
PONIA? PERD-A (2080), Y(200),EZERO, BUIN. THICK. TALFA. APS 
SONMON-MAFO/A(200,59).F(200,1).H(29),T15(200), TSAT.CK, NEL. NSNP, Na 


KN.ICOR(C 200,3) 

LOOPsNEQ-1 

DO 201 Is1,LO0P 

МВ=1+1 

HE=MINOCI+MAXB-1,NEQ) 

DO 281 J=MB.NB 

L=J+2-MB 

ST, LI/ACI.1) 

DO 202 Mm=1.NUEC 
ЕС. ММ) =Р СУ, ММ 2 -БХЕ СТ, ММ) 

MM=MINO (MAXB-L+1,NEQ-J+1 ) 

DO 261 Ks1.mm 

NN=L+K-1 
Хо асл,Козаст.ко-ожаст, мм) 
| DO 2803 I-1,NUEC 
M) F(NEQ, I)-F(NEQ, I2/&(NEQ, 1) 
| DO 894 I-:8.NEQ 

Јемеа-1+1 

К«МІМ2іМЕС-1%1,МАХР) 

DO 204 MM=1,NUEC 
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